function gama_s = ga_s(D,t)
%计算太阳方位角
%参数 %D为距离春分日的天数 t为当地时间（真太阳时）24小时制

global lon %(东)经
global lat %(北)纬
lon=98.5;
lat=39.4;
lat=deg2rad(lat);
sin_sigema=sin(2*pi*D/365)*sin(2*pi*23.45/360);% sigema为太阳赤纬角

alpha_s=alp_s(D,t);
alpha_s=deg2rad(alpha_s);

cos_gama_s=(sin_sigema-(sin(alpha_s)*sin(lat)))/(cos(alpha_s)*cos(lat));

gama_s=acos(cos_gama_s);%计算太阳方位角弧度值
gama_s=real(rad2deg(gama_s));%计算太阳方位角度 取实部

end

